Persistent current and correlation effects in carbon nanotubes 
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The persistent current of interacting electrons in toroidal single- wall carbon nanotubes is evaluated 
within Haldane's concept of topological excitations. The overall pattern of the persistent current 
corresponds to the constant interaction model, whereas the fine structure stems from the electronic 
exchange correlations. 
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The recent breakthrough in the synthesis of a new gen- 
eration ofj-nuantum wires - single-wall carbna nanotubes 
(SWNTs)El and the subsequent observatioro'cl of coherent 
electron transport in this system have initiated a surge nf 
experimental and theoretical activity (see e.g. Refs.cTLl) 
. The investigation of non-Fermi liquid correlation ef- 
fects due to one-dimensional nature of interacting elec- 
trons in SWNTs presents one of the main challenges. The 
signatures of such correlations are often masked by the 
chargiiig effects. Nevertheless, very recent experimental 
resultaJ on the spin structure of the ground state suggest 
the interpretation in terms of electron correlations. 

Generically, carbon nanotubes have linear or curved 
shape. Recently Liu et. alll3 have observed individual 
circular SWNTs and ropes of such nanotubes. The ex- 
perimental data suggests uniform widths of SWNTs and 
does not display the presence of their ends. For this rea- 
son it is plausible that circular SWNTs have the topol- 
ogy of a torus. One can test this conjecture by measur- 
ing the response of circular SWNTs to a weak perpen- 
dicular magnetic fielcL Since SWNTs are typically al- 
most free from defectil, the presence of delocalized elec- 
tronic states in the system should result in a persistent 
current. To provide theoretical support for future ex- 
periments, we analyze the persistent current in toroidal 
SWNTs (TNTs). 

We are aware of two recent paperstiEj where the per- 
sistent current in TNTs has been computed from first 
principles within the single-particle scheme. The results 
are either related to the specific fuUerenes (C576 carhan 
toroidsEJ) , or limited to the special case of half-fillingEj. 
Moreover, both works ignore electron correlations due to 
the Coulomb interaction, whose observable signatures in 
the pattern of the persistent current are investigated in 
this Letter. I_. 

We employ bosonization formalisrr£3 which has proven 
to be effective for analysis ofpflersistent currents in var- 
ious one-dimensional jjBodelaiJ. rFirst, we extend the 
bosonization schemeQ'Ota (see alsoB) to the case of TNTs 
of the "armchair" {N,N) type. The Fermi operators 
^'s for electrons with spin s = ± can be expanded 
near the two crossing (or Fermi) points aK (a = ±, 
K = 47r/3a) of the energy spectrum into right- (d = 
+) and left-moving {d = — ) components, '^s{x) — 
^''°'^^'^asd{x). The pcriodicity of the electronic fields 

s{x + L) ~ '^s{x) results in the following boundary 
condition for the slowly- varying parts, tpasd{x + L) ^ 



4'asd{x)e^^^°'P^^ where p = 0,±1 parametrizes the num- 
ber L/a of elementary cells along the nanotube of length 
L, L/a — in + p (n is an integer). 

Bosonization allows one to express the Fermi oper- 
ators ipasdix) oc e^i't>'>=+d-^''=) tcrms of the bosonic 
fields 9as and (pas obeying the commutation rules 
[0asix),4>a's'{x')] = {ni /2)sign{x - x')Saa' ^ss' ■ The ficlds 
can be decomposed into topological parts and non-zero 
bosonic modes 0as, 4>as- 

9 as (*^ ), (1) 

(bas{x) = + [Jas + 2pa/3)T:x/L + 4> 

as )• (2) 

The pairs of the action-angle operators Jas, 9a] and 
Nasjipas satisfy the canonical commutation relation. 



Nas,(l)] 



(0) 

a' s' 



excitations Nas, Ji 

Ma 



—iSaa'Sss' ■ The topological 
Jas are simply related to the numbers 
i-asd — {Nas + dJas)/2 of cxccss clcctrons at the branch 
a, s, d of the energy spectrum, see inset of Fig. la. Since 
Masd are integers, the sum Nas + Jas must be even (for- 
mally this topological constraint follows from the bound- 
ary condition on t/;-operators). 

Following Refs.LrLl we introduce bosonic fields 9sv{x) 
and (j)sv{x) describing the charge v = + and spin v = — 
excitations in the symmetric (5 = + and antisymmetric 
5 = — modes. 



Osu = [0++ + vO+- + 5{0-+ + vO--)]/2, 



(3) 



where O = 9 {9, 9^°\ N) or (0, (l)'-"\j) and the indices 
in the r.h.s. correspond to a, s. The new fields, 



esAx) = efj + {Ns, + 25s+5,+)t:x/L + 9su{x), (4) 
(l)5u{x) = + iJsi' + (4/3)p55-5i.+)7ra;/L + <f>sAx), (5) 



JSu,9s'u' 



satisfy the commutation relations, [9sw{x),4's'v'(x') 
(7ri/2)sign(x - x')55S'6y, 
—i5ss'6uu'- The topological constraint for Nas + Jas im- 
plies that Yj&v^Sv, Y.5uJsv, Y^v^Sv + Jsv, J2s^Su + 
Jsij, all must be even, whereas J2su ^Su + Jsw = mod 4. 
In addition, the new topological numbers should be either 
all integer or all half-integer. Note that Ntot — 2A^++ + 4, 
2iV++ = 'Y^asd -^^as^is f^c total number of extra elec- 
trons in the systenxj, whereas 2J++ = X^asd '^-^"s'i 
the difference in numbers of right- and left-movers. 
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We. |-|Cpncentrate first on the Luttinger model-like 
termD'cHa of the low-energy Hamiltonian H — Hl + V 
of TNTs, 



(6) 



where the standard interaction parameters Kgu and ve- 
locities vg^ of excitations are introduced, so that Ksv — 1 
for non-interacting electrons. The electro-chemical po- 
tential /i of an electronic reservoir controls the charge 
density p = 2eV0++/7r. The vector potential A of an ex- 
ternal magnetic field is coupled to the persistent current 
/ — — 2e^+_|_/7r. We will neglect the Zeeman term /x^iJ 
which is a factor ^ a/L smaller than the energy scale 
vp/ L of interest. 

The Hamiltonian Hl splits into the bosonic part Hh 
(which has the form (^) with 9 6, (j) ^ cj), and ^ = 
^ = 0) and the topological part. 



2L ^ J^' 



A^i. - 4 ( - - ) 5s+5, 



(7) 



Here = <i>/$o is the magnetic flux <i> through the 
TNT in units of the flux quanta <i>o = 2Tihc/e, and 
= (i4'+4_L/27rw+_|_)/i is a normalized electro-chemical 
potential, whose reference point corresponds to the cross- 
ing of the energy spectrum. The increase of by one cor- 
responds to the addition of an electron to each branch a, 
s, d of the spectrum, so that the properties of the system 
are periodic in with a period of one (the same periodic- 
ity occurs in /$). The Hamiltonian (^ shows additional 
symmetries with respect to changes in sign of the electro- 
chemical potential (/^ ~f^, Masd "> ~Mas^d - 1) or 
the magnetic flux {f^ — > Masd M^as-d)- 

The persistent current / = dF/d^ can be calculated 
by differentiating the free energy F of the system with 
respect to the magnetic flux 



I={ev++K++/L){8U-2{J++)). 



(8) 



Due to the symmetries of the Hamiltonian, the persistent 
current is an even (odd) periodic function of (/*). At 
zero temperature the average ( is determined by the 
ground state, whose map is given in Figs. 1,2. The persis- 
tent current ^ shows saw-tooth dependence on the mag- 
netic flux and changes in a stepwise manner as a function 
of the electro-chemical potential. The amplitude of the 
persistent current is given by /max = 4ew++X++/L w 0.5 
/iA, for ~ vp Ri 8 X 10^ m/s and L ~ 1 /ini. 

This value is by two orders of magnitude larger than tlte 
persistent current measured in GaAs mesoscopic ringsli3. 

The unscreened long-range Coulomb interaction 
strongly influences the forward scattering of electrons 
leading to a large stiffness of the symmetric charge 



modeQS, « 0.2. We will first ignore the sublattice- 
dependent part of the forward scattering as well as the 
backscattering of electrons so thatEl v++ — vf/K++ 
and Ksu = 1, vsi, = vp for the modes {5v) — 
(H — ), ( — h), ( ). Within this approximation, the en- 
ergy spectrum of the topological Hamiltonian (^ is given 
by the sum of the Coulomb and single-particle energies, 
which ccjcresponds to the constant interaction model (see 
e.g. Ref.!). 

The ground state configurations for TNTs withp = 0, 1 
are shown in Fig. 1. Due to the spin degeneracy, the 
eigenstates of the Hamiltonian (^ are characterized by 
four topological numbers. Mad — X)s -^asd, Fig. lb. 
For the nanotubes with p — Q the states of electrons 
moving in the same direction at a = ± are degenerate. 
The system can be described by the two numbers Md = 
^'^ad of extra right- and left-movers, Fig. la. 

Since the Coulomb interaction in SWNTs is strong 
<C 1), the electron number Ntot is determined pri- 
marily by the electro-chemical potential, although in nar- 
row regions it can be controlled by magnetic flux (Fig. 
1). The slope of the ground state borders enables one 
to deduce the value of— the interaction constant 
from experimental datalia. Generically, the changes of 
the ground state with the magnetic flux correspond to 
the jumps of electrons between different branches of the 
spectrum and occur at universal values of magnetic flux 
- = 0,1/2 for p = and = 0,1/6,1/3,1/2 for 
p = ±1. In particular, the jump of an electron at zero flux 
causes the paramagnetic response of TNT. Such param- 
agnetic ground states occur if Ntot ^ 4 mod 8, for p = 0, 
and Ntot 7^ mod 4, for p = ±1 (see Fig. 1). Otherwise, 
the ground state is diamagnetic. 

The sublattice-dependent part of the forward scatter- 
ing and backscattering of electrons lead to the appear- 
ance of an essentially non-Luttinger term V = Vf -\- Vb 
in the Hamiltonian and the renormalization of the pa- 
rameters Ksv, vsi, in Eq. (^). The Luttinger and 
non-Luttinger parts of the Hamiltonian describe intra- 
and interbranch scattering of electrons respectively. The 
derivation of these terrnsJiam a microscopic model has 
been discussed in RefsB'EHla and here we present only 
the results for a generic case away from half-filling. The 
non-Luttinger terms are given by 



Vf^ — ^ dx{cos26i+_cos26i_„ 
27r^a^ J_L/2 

— cos 29 ^ cos 20__ — cos 26 |_ cos 26^ 

Vb = I dx{V{2K) cos 29+- cos 20 



Ti'^a^ J-L/2 
AV{2K) 



(9) 



(10) 



[cos 29 |_ cos 20__ + cos 26 |_ cos 20__ 

-t- COS 26'__ cos 2(/)__]}, 

where the matrix elements V{q) = [V+{q) + V-{q)]/2, 
AV{q) — V+{q) — V-{q) are related to the electron scat- 
tering amplitudes Vp{q) given below, and d ^ 1/K is the 
standard ultraviolet cutoff in the Luttinger model. The 
interaction constants Ksi, = ^/Bs^^JAs^, and velocities of 
excitations vsi^ ~ vpsj As^Bgy (^ are defined byu 
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uation for i? > a givesB AViO) = 0.21, V+{2K) = 0.60, 
V-{2K) = 9.4 X 10-"* in units of ae^/2iTKR (the scatter- 
ing amplitudes AV{Q),V±{2K) increase with decreasing 
the localization radius oq of pa-|Orbital; this phenomeno- 
logical parameter of the modeO is chosen as oq — a/2). 

Perturbation theory with respect to the backscattering 
and the sublattice-dependent p4rt of the forward scatter- 
ing of electrons is applicable i£3 max[Al/(0), V+{2K)] <^ 
2ttvf, see Eqs. dl^), (p^. This condition is equivalent 
to TV = 27ri?/voa ^ c, with c « 0.1 for the param- 
eters listed above, which is safely fulfilled for generic 
SWNTs with TV = 10. The perturbation splits degen- 
erate electronic states n, n' belonging to the same un- 
perturbed energy level i. The splitting occurs already 
in first order and can be estimated from the secular 
equation, det |14,„/ — EiSnn'l = 0. The unperturbed 
states n are characterized by the topological numbers 
Nsi/, Jsv and by the quantum state \...)b of bosonic 
modes. Only the vacuum state |0)(, has to be con- 
sidered at low temperatures T ^ vpjL. The diago- 
nal matrix elements Ynn correspond to the energies of 
the topological excitations (|7|) (we will drop the con- 
stant energy shift due to the renormalization of bosonic 
term (^)). The topological and bosonic parts of non- 
diagonal matrix elements Vnn' (§), (0) can be evalu- 
ated using the relations, e*'^<s>' | Ja'^') = | J^v + l)5ss'Siyi,' , 



{2iT~a/L)^-^, 
we obtain 



\Ns, + l)5ss'5,,', and (Ole^*"^- |0)b = 
22#^.|o)fc = {2T:h/LY/^'-. As a result 



FIG. 1. The ground state of TNTs with p = (a) and p = 1 
(b) in terms of the topological numbers Md = M_) (a) 

and Mad = [M++, M+-, M-+, M ] (b). The ground state 

changes at the sohd lines. We choose = 0.2. Inset: Elec- 
tronic states near the crossing point a = 1 of the spectrum 
for TNTs with p = (circles) —1 (squares) 1 (triangles) at 
/$ = 0. The filling of the electronic states with the energy 
E < 2nvF/3L at the branch a,s,d corresponds to Masd = 0. 
The plot for a — ~1 can be obtained by symmetry with re- 
spect to fc = 0. 
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= 1 



^ 417(0) 


Ay(0) V+{2K) 


(11) 
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Attvf 2'itvf 
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^V+{2K) 
2'kvf ' 


(12) 


^ Ay(o) 




(13) 




2ttvf 



Asy = 1 



Bsu = 1 



Here V+ (g) and V- {q) are the amplitudes of intrasublat- 
tice and intersublattice electron scattering with momen- 
tum transfer (7 = 0, 2K . The forward scattering [q = 0) 
has the strongest amplitude, 1^(0) = 2e^ ln{Rs/ R)/ k, 
where k is an effective dielectric constant of the mfe 
dia (an estimateQ for the parameters of the experiment^ 
gives K — 1.4), i? is the radius of the nanotube, and 
Rs ~ min(L/7r, Z?) characterizes the screening of the 
Coulomb interaction due to a finite length L of the TNT 
and/or the presence of metaUic electrodes at a distance 
DU. The amplitudes AV{0) and V+{2K) decay as 1/R for 
i? > a, whereas V-{2K) < min[Al/(0), V+(2iv:)] due to 
the C3 symmetry of the graphite lattice. Numerical eval- 



{N,J\ / dxcos2es^cos2ds'^'\N',J'] 
L \ Zna ^ 



(14) 



and similar expressions for the other matrix elements. 
Here S... denotes that all topological numbers different 
from Js^, Js'u' should be equal for the initial and final 
states. The last term stems from the integration over x 
in Eqs. (||), (10). It produces an additional constraint 
on the topological numbers, which can be traced back 
to the conservation of momenta of two scattering elec- 
trons. The topological constraint has a somewhat differ- 
ent form, 6pN' +p'j' +2pp',o, for the last term in Eq. 
(|lO|), which contains two field operators for the same 

( ) sector. Let us note that the non-diagonal part 

®)j ( |lO|) of the perturbation does not contain matrix 
elements in the (-I-+) sector. For this reason, perturbed 
ground states are characterized by a well defined topolog- 
ical number J+_|_ (and TV+_|_) which determines the persis- 
tent current (^ at zero temperature. Perturbed ground 
states are shown in Fig. 2. At not very small magnetic 
flux, the perturbation lifts the spin degeneracy of two- 
electron (or two-hole) ground states favoring spin aligned 
configurations (like (2,0) in Fig. 2). With decreasing 
magnetic flux, the "many-particle" ground states (with 
2A^++ ~ 2. ..6 mod 8) experience reconstruction, so that 
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chemical potential. A submicroamp persistent current 
should be observable in a few micrometer long TNTs. 
The Umklapp scattering of electrons on the atomic lat- 
tice (at half- filling), impurities, structural imperfections, 
twiston phonons, etc. may suppress the persistent cur- 
rent and deserves further analysis. 
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FIG. 2. The fine structure of the ground state for TNT 
with p = 0. The numbers M+,AI- of right- and left-movers 
are given in brackets. The parameters AV{0),V±{2K) are 
hsted in the text below Eq. (|l|), K++ = 0.2, and a = 1/K. 
Quantum states in a coherent superposition are denoted by 
double dots. 

both the spin and orbital configurations are changed. 
The reconstruction is observable as a jump of the per- 
sistent current due to the change in numbers of right and 
left movers. The increase of the kinetic energy of new 
ground states is accompanied by the build-up of many- 
electron correlations, which minimizes the total energy. 

For the states (2,0), (2,1), and (2,2) in Fig. 2 the 
electron spins are parallel, which is a signature of the ex- 
change interaction. The non-diagonal terms (||), (|o|) of 
the Hamiltonian do not mix degenerate electron configu- 
rations corresponding to each of these states. Let us note 
that the spin aligned ground states have, been presum- 
ably detected in very recent experimentgj on individual 
linear SWNTs, albeit the data differs substantially from 
the resultsB on ropes of SWNTs. 

The situation is different for the ground states (1,1) 
and (3, 1) (Fig. 2). Each state represents a coherent 
superposition of 4 configurations with antiparallel spins, 
which has the lowest energy due to the interbranch elec- 
tronic exchange allowed by the non-diagonal matrix el- 
ements (|), (PI) of the Hamiltonian. The new ground 
states (1, 1), (272) with even number of electrons are sta- 
ble with respect to a change of sign of the magnetic flux. 
For this reason, TNTs are diamagnetic for even Ntot and 
paramagnetic for odd Ntot, in contrast to the result of 
the constant interaction model. 

In conclusion, we have generalized the bosonization 
formalism for the case of TNTs and evaluated the per- 
sistent current in this system away from half-filling. The 
pattern of the persistent current depends on the num- 
ber of elementary cells along the nanotube modulo 3. 
The overall pattern (Fig. 1) corresponds to the constant 
interaction model, whereas the fine structure (Fig. 2) 
can be explained in terms of electronic exchange cor- 
relations. Even though a system with a fixed electro- 
chemical potential was considered, the results for fixed 
number of particles can be easily obtained from Eq. (|^) 
and Figs. 1, 2 by an appropriate choice of the electro- 
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